Stiffness Exponents for Lattice Spin Glasses in Dimensions d = 3, . . . , 6 
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The stiffness exponents in the glass phase for lattice spin glasses in dimensions d = 3, . . . , 6 are 
determined. To this end, we consider bond-diluted lattices near the T = glass transition point 
p* . This transition for discrete bond distributions occurs just above the bond percolation point pc 
in each dimension. Numerics suggests that both points, Pc and p* , seem to share the same 1/d- 
expansion, at least for several leading orders, each starting with l/(2d). Hence, these lattice graphs 
have average connectivities of a = 2dp > 1 near p* and exact graph-reduction methods become very 
effective in eliminating recursively all spins of connectivity < 3, allowing the treatment of lattices 
of lengths up to L = 30 and with up to 10^ — 10® spins. Using finite-size scaling, data for the 
defect energy width a(AE) over a range of p > p* in each dimension can be combined to reach 
scaling regimes of about one decade in the scaling variable L{p — p*Y . Accordingly, unprecedented 
accuracy is obtained for the stiffness exponents compared to undiluted lattices (p = 1), where scaling 
is far more limited. Surprisingly, scaling corrections typically are more benign for diluted lattices. 
We find in d = 3, . . . , 6 for the stiffness exponents 1/3 — 0.24(1), 1/4 = 0.61(2), 2/5 = 0.88(5), and 
ye = l-l(l)- The result for the upper critical dimension, d„ = 6, suggest a mean-field value of 
ycx> = 1. 

PACS number(s): 05.50.-Hq, 64.60.Cn, 75.10.Nr, 02.60.Pn. 



I. INTRODUCTION 

The stiffness exponent y (often labeled 0) is one of 
the most fundamental quantities to characterize the low- 
temperature state of a disordered spin system . It pro- 
vides an insight into the effect of low-energy excitations 
of such a system 0, ■ A recent study suggested the im- 
portance of this exponent for the scaling corrections of 
many observables in the low-temperature regime Q , and 
it is an essential ingredient to understand the true na- 
ture of the energy landscape of finite-dimensional glasses 

iiQ. 

To illustrate the meaning of the stiffness exponent, one 
my consider an ordinary Ising ferromagnet of size L'^ with 
bonds J = +1, which is well-ordered at T = for rf > 1, 
having periodic boundary conditions. If we make the 
boundary along one spatial direction anti-periodic, the 
system would form an interface of violated bonds be- 
tween mis-aligned spins, which would raise the energy of 
the system by IS.E ^ L'^~^ . This "defect" -energy AE 
provides a measure for the energetic cost of growing a 
domain of overturned spins, which in a ferromagnet sim- 
ply scales with the surface of the domain. In a disordered 
system, say, a spin glass with an equal mix of J = ±1 
couplings, the interface of such a growing domain can 
take advantage of already-frustrated bonds to grow at a 
reduced or even vanishing cost. Defect energies will be 
distributed with zero mean, and the typical range, mea- 
sured by the width of the distribution a{AE), may scale 
like 



a{AE) ~ L» 



(1) 
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Clearly, it must he y < rf— 1, and a bound of ?/ < (d— 1)/2 
has been proposed for spin glass systems generally 2]. 
Particularly, ground states of systems with y < would 
be unstable with respect to spontaneous fluctuations, 
which could grow at no cost, like in the case of the one- 
dimensional ferromagnet where y ~ d — 1 = 0. Such a 
system does not manage to attain an ordered state for 
any finite temperature. Conversely, a positive sign for y 
at T = indicates a finite- temperature transition into an 
ordered regime while its value is a measure of the stabil- 
ity of the ordered state. Furthermore, in a d-dimensional 
family of systems, the marginal value y^^ = provides 
the lower critical dimension dc for such systems. 

Accordingly, there have been many attempts to obtain 
the value of stiffness exponents in finite-dimensional spin 
glasses [lii|ia|nl|ll|ll|lllil|ll|n^, using trans- 
fer matrix, optimization, or renormalization group tech- 
niques. In the early days of spin-glass theory, it was soon 
argued that y < for d < 2 and 2/ > for d > 3 [1 [ill. 
Only recently, though, the stiffness exponent for d = 2, 
below the lower critical dimension, has been improved 
to considerable accuracy, t/2 = —0.282(2) 0,0. There 
has still been little progress in the accurate determination 
of 2/3 in the last 20 years, despite significant increases in 
computational power. It's value is expected to be small 
and positive, and so far has been assumed to be near 
2/3 « 0.19 [Til Il3^. although there have been investiga- 
tions recently pointing to a larger value, such as 0.23 [3 
or 0.27 In some sense, all of these results are con- 
sistent, since they were obtained over exceedingly small 
scaling windows, L — 6, . . . , 12 at the best, and large 
errors have to be assumed . In d = 4 the only value 
reported to date has been 2/4 = 0.64(5) using 

In this paper we use numerical investigations of ±J 
spin glasses on dilute lattices to obtain improved pre- 
dictions for the stiffness exponents in dimensions d — 
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3, . . . , 6. First, we explore such lattices near their bond- 
percolation transition to find a separate transition 
p* > Vr i nto & T — Q spin-glass state, as anticipated by 
Refs. We find that p* becomes ever closer to Pc 

for increasing d, both scaling with l/{2d). Thus, near ei- 
ther transition, bond-diluted lattices have spins with con- 
nectivities distributed near 2dp* « 1. Such sparse graphs 
can be effectively reduced with a set of exact rules that 
eliminate a large fraction of spins, leaving behind a small, 
compact remainder graph that is easier to optimize. The 
increase in the scaling regime with lattice size, in combi- 
nation with finite-size scaling techniques, leads to much 
improved or entirely new predictions for the stiffness ex- 
ponents of low-dimensional lattices. In particular, we 
find that yz = 0.24(1), ^4 = 0.61(2), = 0.88(5), and 
?/6 = 1-1(1) ■ Our value in d = 3 is at the higher end of 
most previous studies and amazingly close to (but dis- 
tinct from) the value obtained with the Migdal-Kadanoff 
approximation, y^^ — 0.25546(3) 17]. The value for 
rf = 4 is consistent with Ref. [lj| and quite below the 
Migdal-Kadanoff value, yf^ = 0.76382(5). The value 
for the upper critical dimension, d„ = 6, seems consistent 
with a recent mean- field prediction of yea = 1, although 
that calculation was based on aspect-ratio scaling ^| . 

In the following section, we discuss the observables that 
our numerical experiments measure, in Sec. IIIII we de- 
scribe the reduction rules for low-connected spins and the 
optimization method use in this study. Sec. IIVI presents 
the results of the experiments for the threshold p*, the 
correlation-length exponent for the glass transition, v* ^ 
and the stiffness exponent y, in each dimension. In 
Sec. El 'we conclude with a discussion regarding the d- 
dependence of y. 



II. DETERMINING STIFFNESS EXPONENTS 

To understand why the accurate determination of these 
stiffness exponents is such a challenging task, it is impor- 
tant to appreciate its complexity: Most numerical studies 
are based on sampling the variance 



(2) 



of the distribution of defect energies /S.E obtained via in- 
verted boundary conditions (or variants thereof p^l. as 
described above. Thus, for an Ising spin glass with pe- 
riodic boundaries, an instance of fixed, random bonds is 
generated, its ground-state energy is determined, then all 
bonds within a hyperplane have their sign reversed and 
the ground-state energy is determined again. The defect 
energy is the often-minute difference between those two 
ground state energies. Then, many such instances of a 
given size L have to be generated to sample the distri- 
bution of and its width (j{/S.E) accurately. Finally, 
a{AE) has to be fitted to Eq. ^ for a sufficiently wide 
range of L in the asymptotic regime. 

The most difficult part of this procedure, limiting the 



range of L that can be achieved, is the accurate deter- 
mination of the ground state energies in the first place. 
While for d < 2 efficient algorithms exist to determine 
ground state energies exactly, and large system sizes can 
be obtained 0, for d > 3 no such algorithm ex- 
ists: Finding ground states is known to be an NP-hard 
optimization problem f23| with the cost of any exact al- 
gorithm likely to rise faster than any power of L. There 
have been a variety of accurate measurements of ground- 
state energies 1211 12^ . |^ using heuristic methods. In 
these measurements small systematic errors in failing to 
obtain a ground state tend to submerge beneath the sta- 
tistical error. In contrast, for the defect energy the exten- 
sive leading-order contributions to the ground states are 
subtracted out, and such systematic failings may surface 
to dominate any statistical errors. Accordingly, system 
sizes that can be approximated with heuristics may turn 
out to be far more limited than one may have anticipated 
based on those previous studies. 

To increase the range of system sizes L without increas- 
ing the optimization problem, we observe that a bond- 
diluted lattice will have the same defect energy scaling 
as a fully connected lattice. Above the finite-size scaling 
window for bond-percolation near pc, the dominant clus- 
ter embedded on the lattice is a compact structure with 
the same long-range properties of the fully connected lat- 
tice. Similarly, the spin-glass problem defined on that 
cluster should exhibit the same long-range behavior as 
the undiluted lattice glass at T = 0, their difference be- 
ing of a short-range geometric nature. Hence, for all bond 
densities p above the scaling window of the T = glass 
transition, Eq. should be applicable. Yet, a spin glass 
on a bond-diluted lattice in turn can be expected to be 
less frustrated, up to the point that frustration fails to 
create long-range correlated behavior. This is certainly 
the case below the bond-percolation transition pc, where 
any defects should remain localized. Thus we focus on 
the regime somewhere above Pc, where the system can 
exhibit spin-glass behavior but where we may take ad- 
vantage of the weakened frustration to optimize larger 
system sizes L. 

As another feature of our new approach, the introduc- 
tion of a new control parameter, the bond density p, per- 
mits a finite size scaling Ansatz. Combining the data for 
all L and p leads to a new variable which has the chance 
of exhibiting scaling over a wider regime than L alone. 
As has been argued in Ref. can make an Ansatz 

of 



a{/^E)^yLyg L{p-p*f 



(3) 



where y ~ 3^o(p~P*)^ refers to the surface tension, which 
must vanish for p ^ p* , and 5 is a scaling function in the 
new scaling variable, x — L [p — p* Y ■ The exponent v* 
describes the divergence of the correlation length for the 
transition into the ordered state at p* . (In the Migdal- 
Kadanoff approximation, it was found that v* is larger 
than V of the percolation transition H^j.) Scale invari- 
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ance at p ^ p* dictates / = yiy*, and in terms of the 
scaling variable x we have 



(4) 



We will use the finite-size scaling relation in Eq. Q to 
analyze our data in Sec. IIV Bl In the following section, 
we describe the new algorithm for spin glasses on dilute 
lattices, which at T = traces out many weakly con- 
nected spins to leave a much reduced remainder graph 
which can be subsequently optimized by other means. 



III. REDUCTION ALGORITHM FOR THE 
ENERGIES 

We will describe the reduction algorithm for spin 
glasses on general sparse graphs at T = in more detail 
elsewhere 24], including its ability to compute the en- 
tropy density and overlap (see also focus here 
exclusively on the reduction rules for the ground state 
energy. We have used these reduction rtdes previously 
for large three-connected Bethe lattices [2^ . These rules 
apply to general Ising spin glass Hamiltonians 



H — ^ ^ '^i^j Xi Xj^ 
<i.J> 



±1) 



(5) 



with any bond distribution P{J), discrete or continuous, 
on arbitrary sparse graphs. Here, we use exclusively a 
± J bond distribution, and bond-diluted hyper-cubic lat- 
tices in d > 3. A Gaussian or any other distribution with 
zero mean and unit variance is expected to yield the same 
value of ?/ jjj . Our preliminary experiments with a Gaus- 
sian distribution have shown faster converging averages 
at a given L, but more persistent scaling corrections for 
large L. 

The reductions effect both spins and bonds, eliminat- 
ing recursively all zero-, one-, two-, and three-connected 
spins and their bonds, but also adding new bonds be- 
tween spins which may or may not have been connected 
previously. These operations eliminate and add terms 
to the expression for the Hamiltonian in Eq. jSJ , leaving 
it form- invariant. Offsets in the energy along the way 
are accounted for by a variable Ho, which is exact for at 
T = 0. 

Rule I: An isolated spin, which does not contribute 
to the sum in Eq. © at all, can be eliminated without 
changing that sum. 

Rule II: A one-connected spin i can be eliminated, 
since its state can always be chosen in accordance with 
its neighboring spin j to satisfy the bond Jij, i. e. in the 
only term in Eq. jS)) relating to Xi, 



XiXj Ji^j ^ l 



(6) 



we can always choose Xi to saturate the bound, which 
is the energetically most favorable state. With that, we 
adjust Ho :— Ho—\Ji.j \ and eliminate the term —Ji.j Xi Xj 
from H. 



Rule HI: A double bond, J-^j and between two 
vertices i and j can be combined to a single bond by 
setting Jij ~ J^^j + or be eliminated entirely, if the 
resulting bond vanishes. This operation is very useful, 
since it lowers the connectivity of i and j at least by 
one. Particular to discrete bond distributions, there is a 
finite probability that the two original bonds cancel each 
other (Jij — 0), which may entirely disconnect i and j 
and reducing their connectivity by two. (Double bonds 
are absent from the original lattice but may arise via the 
recursive application of these reduction rules.) 

Rule IV: For a two-connected spin i, its two terms in 
Eq. lO can be rewritten 

Ji^lXiXi + Ji,2XiX2 = Xi{JisXl -I- Ji,2X2) 

< \JiAXl + Ji,2X2\ (7) 

= Ji,2a;iX2 + AH, 

with 

•^1,2 = -^{\Ji,l + Ji,2\ — \Ji,l — Ji,2\) , 
= i(|J,i + J, 2|-f |J,,i- J,,2|), 

leaving the graph with a new bond J1.2 between spin 1 
and 2, and acquiring an offset Ho := Ho — AH. 

Rule V: A three-connected spin i can be reduced via a 
"star-triangle" relation (see Fig.m): 

Ji,lXiXi + Ji^2XiX2 + Ji^zXiXz) 

= Xi {Ji^ixi + Ji^2X2 + Ji.sxa) 

<\Ji,lXi+ Ji^2X2 + Ji.3X3\ (8) 

= -^i,2a;ia;2 + Ji.as^is^s + -'^2,3a;2a;3 + AH, 

with 

Ji,2 = ~A-B + C + D, Ji,3 = A- B + C - D, 
J2,3 ^ ~A + B + C ~ D, AH = A + B + C + D, 

A — J \ Ji,l — Ji,2 + Ji,3\ 5 ~ i 1*^24 ~ — Ji,3\ j 

= i + ^j,2 + Ji,3\ , D = J \ Ji,l + Ji,2 — Ji,3\ ■ 

The bounds in Eqs. (|6I9|I are saturated for the right 
choice of the spin Xi that links the terms together, thus 
optimizing its alignment with the local field as is required 





1 "2 

FIG. 1: Depiction of the "star-triangle" relation to reduce a 
three-connected spin (xq, center- left). The values for the new 
bonds on the right are obtained in Eq. @. 
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when the remaining graph takes on its ground state. In 
turn, for T > the ehminated spin Xi may not take 
on its own energetically most favorable state to mini- 
mize the free energy of the configuration instead. Hence, 
the reduction algorithm is exact only in determining the 
ground state. 

Reducing four- and higher-connected spins would lead 
to new bonds that connect more than 2 spins, creating in 
general a hyper-graph with multi-spin interaction terms. 
For instance, a term in H connecting a spin ao to four 
other spins would be replace by one term connecting all 
four, six terms mutually connecting the four neighbors in 
all possible pairs, and an energy offset[32|- While such a 
strategy may be useful, we will confine ourselves here to 
reductions producing only new two-spin interactions. 

It is important that these rules are applied recursively 
and in the given order. That is, one may only apply 
Rule II after there are no more spins reducible by Rule 
/, apply Rule III only after both. Rule I and Rule II, 
have been exhausted, etc. And after the application of 
any higher rule, it needs to be checked if structures have 
been generated to which any lower rule may now apply. 
For example, the recursion may have generated a spin 
that is two-connected, but via a double bond to a single 
other spin. Applying Rule IVto that spin before Rule III 
would lead to the other spin having a bond onto itself, 
a problematic situation for which we have no rule. In 
any event, even if we had provided more rules for all 
eventualities, it is still far more efficient to first reduce 
the lowest connected spin at any one time. 

After all these rules have been exhausted, the origi- 
nal lattice graph is either completely reduced (which is 
almost certainly the case for p < pc), in which case Ho 
provides the exact ground state energy already, or we are 
left with a much reduced, compact graph in which no spin 
has less than four connections. Note that bonds in the 
remainder graph may have properties uncharacteristic of 
the original bond distribution. For example, ± J-bonds 
may have combined to bonds of any integer multiple of 
J (e. g. via Rule III). Here, we obtain the ground state 
energy of the reduced graph with the extremal optimiza- 
tion heuristic 23], which together with Hq provides a 
very accurate approximation to the ground state energy 
of the original diluted lattice instance. Clearly, we could 
have just as well used other heuristics or exact methods 
to treat the remainder graph. 



IV. NUMERICAL RESULTS 

The following data was obtained during a window of 
about two months on a cluster of 15 Pentium4 PC run- 
ning at 2.4GHz with 256MB of RAM. 

The runtime of the EO heuristic was fixed to grow as 
with the number n of spin variables in the re- 
mainder graph after the reduction had been applied. In 
Ref. [2^ it was found that typically 0{n'^) updates for in- 
stances up to n ~ 10'^ were needed to obtain consistently 



TABLE I: List of the bond-density tliresholds on liyper-cubic 
lattices for percolation pc (tafcen from Ref. |2^ ') and for the 
T = 0-transition into a spin glass state, p* , as determined 
from Figs. |21 The values of p* are dependent on the bond 
distribution which is ±J here. 



d 


Pc P* 


3 


0.2488 0.272(1) 


4 


0.160130 0.1655(5) 


5 


0.118174 0.1204(2) 


6 


0.0942 0.0952(2) 



reproducible ground state energies. Since we are aiming 
at much larger statistics and typically smaller instances 
in the present study, we opted for a more limited run- 
time. Instead, an adaptive multiple restart system was 
used such that for each instance at least 3 runs from fresh 
random initial spin configurations were undertaken. If a 
new best-so-far configuration is found in run r, at least 
a total of 2r restarts would be applied to these more 
demanding instances |2^. For instances with n > 700 
apparent inaccuracies in sampling the difference between 
ground state energies, AE, become noticeable. 

For highly connected graphs with few spins to reduce, 
local search with the EO heuristic dominated by far the 
computational time. Our implementation of the reduc- 
tion algorithm, originally conceived with d = 3 lattices 
with up to L = 30 in mind, started contributing sig- 
nificantly to the computational cost for instances with 
^ 20^, hence most noticeably in the study of p* in 
d = 5 and 6. 



A. Determination of p* 

In Ref. |0| it was shown that spin glasses on diluted 
lattices may possess a distinct critical point p* in their 
bond fraction, which arises from the (purely topological) 
percolation threshold Pc of the lattice in conjunction with 
a discrete distribution of the bond weights P{J)- Clearly, 
no long-range correlated state can arise below pc- A crit- 
ical point distinct from percolation, p* > pc, emerges 
when such an ordered state above Pc remains suppressed 
due to collaborative effects between bonds [3| (see Rule 
III in Sec. IIIip . Just above Pc, the infinite bond-cluster 
is very filamentary and may easily be decomposed into 
finite components through such collaborative effects, in- 
volving a small number of bonds along narrow "bridges" 
between those components. Thus, to observe the onset of 
glassy properties on a dilute lattice, we have to cross an- 
other threshold p* > Pc first. In Ref. we were able to 
locate p* for the Migdal-Kadanoff lattice in accordance 
with theory [l^ by using the defect energy scaling from 
Eq. : For all p > p* the stiffness exponent y eventu- 
ally took on its p = 1 value, while for any p < p* defect 
energies diminished rapidly for increasing L. 

In each dimension, we have run the above algorithm on 
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FIG. 2: Plot on a logarithmic scale of the variance a{AE) of 
the defect energy as a function of systems size L for various 
bond fractions p > in d = 3 to 6. In each case, a{AE) drops 
to zero rapidly for increasing L at smaller p, but turns around 
and rises for larger p, indicative of a nontrivial glassy state 
at low T. Near p* , a{AE) undergoes ever longer transients. 
The values for the thresholds p* as suggested by the plots are 
Usted in Tab. |T| 
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FIG. 3: Plot of the known bond-percolation thresholds pc, 
the T = glass transition thresholds p* determined in Fig. |21 
and the 1 /d-expansion of pc in Eq. JSJ as a function of d. The 
insert shows ln(p* —pc)/ln{2d) vs. l/ln(2ci), extrapolating 
seemingly toward a > 4 as in Eq. (Unj. 



a large number of graphs (about 10^ — 10^ for each L and 
p) for p increasing from pc in small steps. For each given 
p, L increased until it seemed clear that (7{AE) would 
either drop or rise for good. In this way, we bracket- 
in p*, as shown in Figs. [3 Both, the bond-percolation 
thresholds Pc, taken from Ref. [23|, and our results for p* 
are listed in Tab. |I| 

It is interesting to compare the values of p* to those of 
Pc for increasing dimension d. In Fig.|31we plot both, pc 
and p* , as a function of d, together with the prediction 
of the three-term 1/d expansion [2^ 



1 

2d 



Pc 

The difference p* 
suming 

P -Pc 



(9) 



Pc clearly decreases for d oo. As- 



(2d)-" {d~^oo), 



(10) 



we plotted ln{p* —pc)/\n{2d) vs. l/ln(2fi) in the in- 
sert of Fig. O to extrapolate for a. This crude extrap- 
olation suggests a > 4, so that p* may share the 1/d- 
expansion of Pc in Eq. 0, at least up to the given order. 
In any case, a bond-diluted lattice system with discrete 
±J bonds enters its spin glass phase at an average con- 
nectivity 2dp* ~ 1, and the reduction methods outlined 
in Sec. lIIII should be very effective in any sufficiently large 
dimension for p ^ p* ■ 

The value of p* is distribution-dependent i ^^'^ the 
values determined here and listed in Tab. ^ result from 
discrete ±J-bonds. It is expected that p* — pc for any 
continuous distribution. The precise values for p* , while 
interesting in their own right, are not important for the 
following discussion of the defect energy scaling. We 
merely need to ensure a selection of bond densities suffi- 
ciently above p* , where we would expect Eq. ^ to hold, 
and sufficiently close to p* for an effective application of 
the reduction rules in Sec. IIIII 
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B. Determination of Defect Energy Scaling 

We have conducted extensive numerical experiments 
to extract the asymptotic scaUng of a{AE) for a many 
conveniently chosen bond densities p, especially in d = 3, 
but also in higher dimensions, up to the upper critical 
dimension ci = 6 Q]. As mentioned in Sec.^ an appro- 
priate choice of p is crucial to ensure a good compromise 
between maximal algorithmic performance (for smaller 
p > p*) and minimal scaling corrections (for larger p) 
that maximizes the actual scaling window. While we can 
estimate the effect of p on the performance of our al- 
gorithm, we have a-priori no information about scaling 
corrections. We will see that scaling corrections are in- 
deed large for p — > p*. Yet, as luck will have it, they 
diminish rapidly for intermediate values of p and again 
increase for p — > 1 (at least in lower dimensions, where 
this limit was considered). 

For the study of p* , in principle very large system sizes 
can be reached due to the complete reduction of very 
sparse graphs. Since optimizing the spin glass on the re- 
mainder graph is an NP-hard problem, we have obtained 
more limited maximal system sizes above p* , dependent 
of the bond density p. We obtained sizes ranging up to 
i = 30 at p = 0.28 to i = 9 at p = 1 in d = 3, L = 15 at 
p = 0.18 to L = 5 at p = 1 for d = 4, L = 13 at p = 0.125 
to L — 5 at p = 0.22 in d = 5, and L — 9 a,t p = 0.1 to 
L — A at p = 0.17 in d = 6. For each choice of L and 
p, we have sampled the defect energy distribution with 
at least N > 10^ instances, then determined its variance 
a{AE). For each data point for a{AE) we estimated its 
error bar as 7/^/N. In Figs. 01 we plot all the data for 
each dimension simply according to Eq. JQ), on a loga- 
rithmic scale. For most sets of graphs, a scaling regime 
(linear on this scale) is visible. Yet, various deviations 
from scaling can be observed. Clearly, each sequence of 
points should exhibit some form of finite size corrections 
to scaling for smaller L. For large L, the inability to 
determine defect energies correctly (according to the dis- 
cussion in Sec. will inevitably lead to a systematic 
bias in a. Some data sets did not exhibit any discernible 
scaling regime whatsoever, most notably our data set for 
the undiluted lattice in d = 3. 

To obtain an optimal scaling collapse of the data, we 
focus on the data inside the scaling regime for each set. 
To this end, we chose for each data set a lower cut in L by 
inspection. An appropriate high-end cut is introduced by 
eliminating all data points for which the remainder graph 
had a typical size of > 700 spins; at that point the EO 
heuristic (within the supplied runtime) seems to fail in 
determining defect energies with sufficient accuracy. All 
the remaining data points for L and p are fitted to a 
four-parameter scaling form. 



(11) 




10 

L 




i. e. approximating the scaling function g{x) from Eq. Q 
merely by unity, its leading behavior for large argument. 



FIG. 4: Plot on a logarithmic scale of the width a of the defect 
energy distribution as a function of system size L. From top 
to bottom, the data for dimensions d = 3, d = 4, d = 5, and 
d = 6 is shown. The data is grouped into sets (connected by 
lines) parameterized by the bond density p. Most sets show 
a distinct scaling regime as indicated by Eq. Q for a rang of 
L above finite scaling corrections but below failing accuracy 
in the optimization heuristic. 
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TABLE 11: List of the fitted values for the critical bond- 
density p*,the correlation-length exponent v* , the surface 
tension ^o, and the stiffness exponent y in each dimension 
3 < d < 6. Included is also the Q- value for each fit. The 
values for p* here are bound to be less accurate than those 
directly determined in Sec. EVU but are consistent. In con- 
trast, the values for y are quite stable. 



P 



yo y 



Q (DoF) 



0.2706 f.l7 2.37 0.239 f.OO (92) 

0.1699 0.60 2.43 0.610 0.00 (47) 

0.1217 0.50 3.05 0.876 0.86 (48) 

0.0959 0.44 3.87 1.103 0.02 (46) 



Unfortunately, we have no knowledge of the functional 
form for finite-size corrections, making the low-L cut on 
the data a necessity. The fitted values for this and the 
other fitting constants (p*, v* , and y) are hsted in Tab. ITU 
Holding p* fixed at the independently determined values 
from Tab. ^reduces the variance in the remaining fitting 
parameters without changing much in their quoted val- 
ues. Using the parameters of this fit, we re-plot only the 
data from the scaling regime in each dimension in Figs.|Sl 

In each case, a convincing scaling collapse is obtained. 
Clearly, our data for d = 3 is not only the most extensive, 
but also happens to scale over nearly two decades with- 
out any discernible deviation or trend away from pure 
power-law scaling that would betray any systematic bias 
or lack of asymptotic behavior. This justifies a certain de- 
gree of confidence to project — 0.24(1) for the scaling 
exponent where the quoted error is based on the uncer- 
tainty in the fit. Troublesome is the observation that the 
data for the undiluted lattice (p = 1) never reaches the 
scaling regime (see Figs. |5| top). This may be in accor- 
dance with the observation of Ref. , which found very 
long transients in a similar study on undiluted Migdal- 
Kadanoff lattices (see also Ref. [13 )j or similar findings 
for undiluted lattices [s^l- In our data, systematic er- 
rors in sampling ground states seem to set in for large L, 
before any scaling regime is reached at all. 

For increasing dimension d, accessible scaling regimes 
become shorter, leading to more difficulty in determin- 
ing an accurate fit of the power law. In d = 4 we can 
still claim scaling for about a decade in the scaling vari- 
able, justifying a prediction of y^ = 0.61(2). In d = 5 
and d = 6, we only reach scaling windows significantly 
shorter than a decade. Luckily, yd increases with increas- 
ing d, thus larger absolute errors still result in acceptable 
relative errors, and we predict from the fits in Figs.|31that 
2/5 = 0.88(5) and y^ = 1.1(1). 



V. CONCLUSIONS 

We have used the combined effort of an exact reduction 
method and an efficient heuristic to determine the defect 
energy distribution for ± J-spin glasses on bond-diluted 




L(P-P y 




L(p-p y 

FIG. 5: Scaling plot of the data from Figs. |1] for a, fitted to 
Eq. im as a function of the scaling variable x = L(p — p*Y . 
Data above or below the scaling regime in each set from Figs.|l| 
was cut. From top to bottom, the scaling collapse of the data 
for dimensions d = 3 to 6 is shown. The lines represent a 
power-law fit of the collapsed data which provides an accurate 
determination of the stiffness exponent y in each dimension. 
For d = 3 (top), we have also included the data for p = 1, 
which does not appear to connect to the scaling regime. 
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lattices in low dimensions. A subsequent finite size scal- 
ing fit of the data allowed us to extract the stiffness expo- 
nents in these dimensions to within 4% to 10% accuracy. 
Our approach also allowed the determination of a variety 
of other observables associated with the T — transition 
into a glassy state at a bond-density p* for d > 3. We 
hope that the methods introduced here may be applicable 
as well to the treatment of other open questions regarding 
the low-temperature state of spin glass systems |j, 0] . 

Our value of = 0.24(1) in d = 3 is near the higher 
end of previous estimates varying between 0.19 and 0.27, 
while ?/4 = 0.61(2) in d = 4 overlaps with a previous 
result of 0.64(5) from Ref. U]. There has been no pre- 
vious value for d = 5, but for d = 6, the upper critical 
dimension, our value of j/g — 1-1(1) suggest a mean field 
result of Uoo = 1- There has been no previous determi- 



nation of the exponents v*, except that it is bound to 
exceed the value of v for percolation and that its 
mean field value for d > 6 should he = 1/2. In light 
of the fact that i' 0.875, 0.68, 0.57, and 0.5 for bond 
percolation in d = 3 to 6 [31] . most of our fitted values 
for i^* do not seem to satisfy these expectations, which is 
easily explained with their poor accuracy. For instance, 
i/* > 1/2 should hold, so it appears that the fitted values 
of ly* are generally too low. 

I would like to thank A. Percus, R. Palmer, M. 
Palassini, and A. Hartmann for helpful discussions, and 
our IT staff for providing access to our student comput- 
ing lab. This project is supported by grant 0312510 from 
the Division of Materials Research at the National Sci- 
ence Foundation. 
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